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The magnetic bipolar transistor (MBT) is a bipolar junction transistor with an equilibrium and 
nonequilibrium spin (magnetization) in the emitter, base, or collector. The low-injection theory 
of spin-polarized transport through MBTs and of a more general case of an array of magnetic p-n 
junctions is developed and illustrated on several important cases. Two main physical phenomena are 
discussed: electrical spin injection and spin control of current amplification (magnetoamplification) . 
It is shown that a source spin can be injected from the emitter to the collector. If the base of an MBT 
has an equilibrium magnetization, the spin can be injected from the base to the collector by intrinsic 
spin injection. The resulting spin accumulation in the collector is proportional to exp(qVbefkBT), 
where q is the proton charge, Vb e is the bias in the emitter-base junction, and ksT is the thermal 
energy. To control the electrical current through MBTs both the equilibrium and the nonequilibrium 
spin can be employed. The equilibrium spin controls the magnitude of the equilibrium electron and 
hole densities, thereby controlling the currents. Increasing the equilibrium spin polarization of the 
base (emitter) increases (decreases) the current amplification. If there is a nonequilibrium spin in 
the emitter, and the base or the emitter has an equilibrium spin, a spin-valve effect can lead to a 
giant magnetoamplification effect, where the current amplifications for the parallel and antiparallel 
orientations of the the equilibrium and nonequilibrium spins differ significantly. The theory is 
elucidated using qualitative analyses and is illustrated on an MBT example with generic materials 
parameters. 

PACS numbers: 72.25.Dc,72.25.Mk 
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I. INTRODUCTION 

Integrating charge and spin properties of semiconduc- 
tors is the central goal of semiconductor spintronics 1 
whose prospect has been fueled by the experimen- 
tal demonstration of electrical spin injection into 
semiconductorS ) 2 i 3 i 4 i 5 : 6 as well as by the discovery of 
III-V ferromagnetic semiconductors^ (Eu-based ferro- 
magnetic semiconductors have even been used earlier as 
effective spin fllter a 9 i 10 ) and observations of relatively 
long spin relaxation times^i*^ Many important advances 
have already been made toward an efficient spin control 
of electrical current in semiconductors, and, vice versa, 
control of magnetism by electrical means. Recent ex- 
amples include a control of ferromagnetism by incident 
lighfci&iai or by gate voltageji£*i& spin injection induced 
magnetoresistance in nonmagnetic semiconductors^ or 
the spin-galvanic effect^ 

Transistors are naturally suited for spin control of elec- 
trical currents since the three regions, emitter, base, and 
collector, can serve as a spin injector, transport medium, 
and spin detector, respectively. There has been remark- 
able experimental progress with employing hybrid fer- 
romagnetic metal and semiconductor structures as the 
hot-electron transistors^ the magnetoresistance of such 



5 MBTs were first proposed in Ref. |26| (see also 
where iwe analyzed spin injection and current 



transistors can be as large as 3400%-° and they can 
be used as effective spin injectors^. The theoretical 
proposals for spin transistors focus largely on the field- 
effect systems In this paper we analyze magnetic 
bipolar transistors (MBTs) which are conventional (spin- 
unpolarized) bipolar junction transistors^ with added 
spin 
Ref. |2_ 

amplification of npn MBTs with a source spin^ Spe- 
cial cases of MBTs without a source spin were recently 
studied by Flatte et al^£ who calculated the spin cur- 
rent polarization in a magnetic-base npn MBT, and by 
Lebedeva and KuivalainenSi, who calculated the cur- 
rent amplification in a magnetic emitter pnp MBT. Of 
spin transistors, the closest one to MBT is the so called 
SPICE (spin polarized injection current emitter) ^ Q which 
employs ferromagnetic metals in the emitter and base- 
collector regions. One of the principle drives for propos- 
ing all-semiconductor spin transistors is the possibility of 
controlling current amplification by spin. 

MBTs integrate ferro (magnetic) and nonmagnetic 
semiconductors in the usual bipolar junction transis- 
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hybrid ferromagnet/semiconductor heterostructures are 
currently an active area of research^ The potential 
of ferromagnetic semiconductors for bipolar devices has 
been shown already in Ref. |3j| where a ferromagnetic 
diode was presented^ More recently (Ga,Mn)As/GaAs 
p-n heterojunctions have been fabricated^ and electri- 
cal spin injection through magnetic bipolar tunnel junc- 



tions has been demonstrated 35,36 showing up to w 80% 
injected electron density spin polarization at 4.6 K— Fi- 
nally, in Ref. a CoMn doped p-Ge and an n-Ge were 
put together to form a ferromagnetic p-n heterojunction 
which showed magnetization dependent current rectifi- 
cation, with up to 97% electrical current variations due 
to the applied magnetic field. Such hybrid junctions can 
also be used for MBTs, where the requirement is that the 
magnetic region has a sizeable equilibrium spin polariza- 
tion (say, 10%). This polarization can be provided by 
the exchange splitting in ferromagnetic semiconductors, 
or by the large Zeeman splitting of dilute magnetic semi- 
conductors. For example, Zeeman splitting can be signif- 
icantly enhanced by large effective g-factors in magneti- 
cally doped (\g\ ~ 500 in Cdo.95Mno.05Se at low tempera- 
tures) or in narrow band gap semiconductors (\g\ ~ 50 in 
InSVj^ even at room temperature). Another possibility 
would be to use a ferromagnetic semiconductor slightly 
above its Curie temperature, a regime also expected to 
give large g-factors. However, before there is an addi- 
tional progress in fabricating junctions using reported 
room temperature ferromagnetic semiconductors (for ex- 
ample, (Zn,Cr)Ta^), the demonstration of the operation 
of MBTs will likely be limited to temperatures below 
~ 150 K. 58 Room temperature MBT is certainly a great 
challenge. 

We formulate here a fully analytical theory of spin- 
polarized transport through MBTs in the small bias (low 
injection) regime, where the injected carrier densities are 
smaller than the equilibrium ones. The theory uses the 
generalized Shockley model for the spin-polarized trans- 
port through magnetic p-n junctions^ as well as the the- 
ory of conventional bipolar junction transistors, as devel- 
oped by ShockleyiS 5 ^ 3 ^ Our theory can thus be viewed as 
a generalized Shockley theory of bipolar transistors. Two 
novel phenomena are studied in detail: electrical spin in- 
jection from the emitter to the collector, and spin control 
of the current amplification (also called gain). Electrical 
spin injection is shown to be effective in the amplifica- 
tion mode of the transistor, the mode where the tran- 
sistor amplifies current. Spin control can be achieved 
by modifying both the equilibrium and nonequilibrium 
spin, since both can modify the electrical current. The 
control by the equilibrium spin (what we call the magne- 
toamplification effect) results from the dependence of the 
equilibrium minority carrier density on the equilibrium 
spin polarization, while the control by the nonequilib- 
rium spin (what we call the giant magnetoamplification 
effect) controls the current via the spin-charge coupling 
of the Silsbee- Johnson typeM^ 

We first describe the model of MBTs in Sec. [H] and 
formulate the analytical theory in Sec. IIIII leaving the 
formal aspects of the theory for the Appendix. We then 
apply the theory to study electrical spin injection through 
MBTs in Sec. IIVI and spin control of the current ampli- 
fication in Sec.[V] where we also discuss the spin current 
in MBTs. 
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FIG. 1: Scheme of a magnetic npn transistor in equilibrium. 
The conduction band E c is populated mostly with electrons 
(filled circles) in the emitter and collector. In this example the 
base is magnetic with the conduction band spin split by 2qC,b- 
The valence band E v , separated by the energy gap E g from 
the conduction band is populated mostly by holes (empty cir- 
cles) in the base. The small population of electrons in the base 
has an equilibrium spin polarization aob = tanh(g("i,/fcsr), 
holes are assumed spin unpolarized. The electron spin is in- 
dicated by the dark (spin up) and light (spin down) shadings. 
The Fermi level (chemical potential) Ef is uniform. Between 
the bulk regions, built-in potentials V be and V b h c are formed 
defining the depletion layers (shaded) in the base-emitter (6- 
e) and base collector (b-c) junctions, respectively. Finally, 
labels "source", e, be, be, and c stand for the regions at which 
they are shown. For example, be is the region in the base at 
the boundary with the depletion layer. 



II. MODEL 

A conventional, spin-unpolarized bipolar npn 
transistor— consists of n, p, and n regions connected in 
series (consult Figs. ^ and [2J. Typically the n region 
with the higher donor doping is called the emitter, the 
one with the lower doping the collector. The base is the 
p region (doped with acceptors) sandwiched in between. 
The most useful mode of operation of the transistor is 
the so called amplification (also forward active) mode, 
where the emitter-base (b-e) junction is forward biased, 
so that the electrons are easily injected into the base. 
Together with the opposite flow of holes, they form the 
emitter current j e . The electrons injected into the base 
diffuse towards the collector. The base-collector (b-c) 
junction is reverse biased. This means that any electron 
reaching the junction from the base is swept by the elec- 
tric field in the depletion layer to the collector, forming 
the electron current (holes's contribution is negligible). 
The base current is the difference jb — j c — je- This 
difference comes from two sources. First, from the hole 
current which is present in the emitter but not in the 
collector. Second, from the electron-hole recombination 
in the base which diminishes the number of electrons 
that make it to the collector. These two factors form 
the generally small jf,. The current amplification (3 is 
defined as the ratio of the large collector current to 
the small base current. For practical transistors /3 is of 
order 100, meaning that small variations in jb (input 
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FIG. 2: Scheme of the magnetic npn transistor from Fig.Qin 
the amplification mode. The b-e junction is forward biased 
with Vbe > 0, lowering the barrier and reducing the depletion 
layer width. The b-c junction is reverse biased with Vbc < 0, 
raising the barrier and increasing the depletion layer width. 
The corresponding changes to the Fermi level Ef are indi- 
cated. The emitter has a spin source, indicated here by the 
incident circularly polarized light generating nonequilibrium 
electron spin well within the spin diffusion length L 3 from 
the b-e depletion layer. The electron and hole flow gives the 
emitter (j e ), base (jb), and collector (j c ) charge currents. The 
electron-hole recombination is depicted by the dashed lines. 
Also shown are the effective widths of the emitter (w e ), base 
(wb), and collector (w c ). 



signal) lead to large variations in j c (output signal). To 
maximize the gain one needs to (i) minimize the relative 
contribution of holes in j e , or (ii) inhibit the electron- 
hole recombination in the base. Typically silicon is used 
to make bipolar transistors, since the indirect gap makes 
it a poor material for the electron-hole recombination. 
We will show below that MBTs allow spin control of the 
gain by realizing (i). 

The magnetic bipolar transistor is a bipolar junction 
transistor with equilibrium spin due to spin-split carrier 
bands, as well as with a nonequilibrium source spin intro- 
duced, for example, by external electrical spin injection 
or optical orientation.-^ The equilibrium spin can be a re- 
sult of the Zeeman splitting in an applied magnetic field 
or of the exchange splitting due to ferromagnetic semi- 
conductors integrated into the device structure. For our 
purposes the equilibrium spin splitting should be on the 
order of thermal energy for the spin-charge coupling dis- 
cussed below to be significant. If no equilibrium spin is 
present, this restriction becomes irrelevant, but the spin 
effects are limited to electrical spin injection. 

An MBT in equilibrium is described in Fig. ^ The 
base is magnetic, with the spin splitting 2q£b. In the 
emitter the majority carriers are electrons whose num- 
ber is essentially Nd e , the donor density. Similarly in the 
collector, where the donor density is Nd c - Holes are the 
minority carriers in the two regions. The base is doped 
with N a b acceptors. Holes (electrons) are the majority 
(minority) carriers. We assume that only electrons are 
spin polarized. The inclusion of the hole spin polariza- 



tion is straightforward and adds no new physics to our 
considerations. Furthermore, in many important semi- 
conductors (such as GaAs) holes lose their spin orienta- 
tion very fast^ and indeed can be treated as unpolarized. 
Note that the electron density is n = rij + nj_, the (elec- 
tron) spin density is s — n-\— n±, and the spin polarization 
is a — s/n. (Spin polarization of the charge current is 
studied in Sec. ESI ) 

The equilibrium in MBTs can be disturbed by apply- 
ing a bias as well as by introducing a nonequilibrium 
source spin. Figure [21 depicts the nonequilibrium physics 
and introduces the relevant notation. We assume that 
the source spin is injected into the emitter within the 
spin diffusion length from the b-e depletion layer so that 
enough spin can diffuse to the base. At the b-e depletion 
layer the electrons feel a spin-dependent barrier: in Fig. [21 
the barrier is smaller (larger) for the spin up (down) elec- 
trons. As in the conventional bipolar transistors there is 
a significant accumulation of the minority carriers around 
the forward biased b-e depletion layer, while there are few 
carriers around the reverse biased b-c layer. The widths 
w of the bulk regions depend on the applied voltages as 
well as on the equilibrium spin polarization^! 

We assume that the electron-hole recombination, oc- 
curring mostly in the emitter and the base is spin in- 
dependent, a reasonable approximation for unpolarized 
holes. We also assume that the spin splitting is uniform 
in the bulk regions, eliminating magnetic drift (magnetic 
drift in semiconductors is discussed in Ref. 41). Our 
other assumptions are those of the standard Shockley 
theoryi2£i2i: Temperature is large enough (say, T > 50K ) 
for all the donors and acceptors to be ionized; The car- 
riers obey the nondegenerate Boltzmann statistics; The 
injected minority carrier densities are much smaller than 
the equilibrium densities; The electric fields in the bulk 
regions are small eliminating electrical drift. Further- 
more, we neglect the carrier recombination and spin re- 
laxation inside the depletion layers. These effects are 
important at very low biases and are not relevant for our 
observations, although may by themselves lead to nice 
physics. Finally, the contacts with the external electrodes 
are ohmic, maintaining the carrier (but not necessarily 
spin) densities in equilibrium. 



III. THEORY 

We generalize the Shockley theory of bipolar transis- 
tors to include spin. The theory is valid in the small bias 
regime and is applicable to any operational mode of the 
transistor, not only to the amplification regime. Physi- 
cally, the theory describes electron and hole carrier and 
spin diffusion in the bulk regions, limited by the electron- 
hole recombination and spin relaxation. The depletion 
layers provide only boundary conditions for the diffusion, 
by connecting the charge and spin currents in the adja- 
cent regions. The most essential assumption is that the 
spin-resolved chemical potentials remain constant across 
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FIG. 3: Magnetic p-n junction. The input densities are the 
carrier and spin densities at the end of the p region: n v and 
s p , as well as at the end of the ti region* Sn- The densities to 
be calculated are at the edge of the the depletion layer: sl in 
the p and sr in the n region. Also indicated are the effective 
widths of the two bulk regions. 



the depletion layers. 

The transistor is viewed as two p-n junctions in se- 
ries. The minority carrier density in each junction (b-e 
and b-c) is determined by the bias voltage across each 
junction. In MBTs the densities are determined also by 
the spin polarization, which needs to be calculated self- 
consistently, as is explained below. Within the limits of 
the theory it is enough to know the minority electron 
densities n& e and nb c to determine the electron charge 
currents, and p e and p c to determine the hole charge cur- 
rents (see Fig. for labeling the regions). We divide the 
presentation of the theory into two steps. First, we re- 
call the main results of the generalized Shockley theory 
of magnetic p-n junctions^ and second, we use these 
results to formulate a theory of a series of magnetic p- 
n junctions and solve it for npn MBT. The first step is 
necessary to also understand our qualitative analyses of 
the transistor operations in the amplification mode. The 
second step, which is rather technical, is left for the Ap- 
pendix. 

In the following we present selected results of the the- 
ory of magnetic p-n junctions. The notation, which is 
easily adapted for use in MBTs, is summarized in Fig. [3] 
Both the p and n regions are in general magnetic, biased 
with voltage V . The nonohmic contact (to simulate the 
conditions at the base of a transistor) at the p region 
maintains nonequilibrium electron n p and spin s p den- 
sities. Similarly, there are nonequilibrium densities til 
and sl at the left of the depletion layer. In the n region 
electrons are the majority carriers whose densities can 
be considered fixed by the donor density jVjj. However, 
the spin density can vary, being sr at the right of the 
depletion layer and s n at the contact with the external 
ohmic electrode. We use subscript to denote equilib- 
rium quantities. The equilibrium minority densities are 
no p (electrons in p), po n (holes in n), and sq p (electron 
spin in p) . These densities are uniform across the corre- 



sponding regions. The equilibrium density spin polariza- 
tions in the n and p regions are «on and ao p , respectively. 
These are also uniform. We denote the excess (injected) 
quantities by 5. For example, 6s l — sl — sq p . We next 
denote by L np the electron diffusion length in p, and by 
L sn and L sp the electron spin diffusion length in n and 
p, respectively. Finally, D np (D nn ) stand for the electron 
diffusion coefficients in n (p). Similarly for holes. 

The spin injection efficiency in magnetic p-n junctions 
is measured by an — Ssji/Nd, where^i 

Ssr = jo6s n + ~fi(5s p - a 0p dn p ) + ~/ 2 aoL$n p (1) 
- 72 cosh (w p /L np ) sol (e« y / feBT - l) ; (2) 

the transport/geometry 7 factors are 



1 



To 
7i 
72 



cosh(w„/L s „) ' 
D np L sn \ tanh(ui„/L s „) 



DnnL sp J smh(w p /L sp ) 
D np L sn \ tanh(w„/L s „) 



D nn L np J smh(w p /L 

np J 



(3) 

(4) 
(5) 



Equation JIJ is accurate up to the terms of the relative 
order of noexp(qV/kBT)/Nd- While such terms can be 
safely neglected when dealing with the spin and carrier 
densities, they must be included when calculating the 
spin current in the n region, where a difference between 
two small spin densities of the same order needs to be 
evaluated. (These terms are not presented in Ref. l4ll ) 
The exact formula for the injected spin density 5sr can 
be cast in the form of Eq. , but with the coefficients 
7 divided by the factor 1 + v: 



7 -> 7/(1 + v), 



where 



n 0p 



e qV/k B T 



73«0p 



— [71 cosh (w p /L sp ) — 

«0p — C*0n , 



*0p 



*0n 



1 



oc, 



(6) 

(7) 
(8) 



On 



Typically v is a number smaller than 0.1, so the correc- 
tions to the spin injected density are upmost 10%. Know- 
ing San we can calculate the injected minority densities 
Sul and 5s l: 



Sn L = n Qp 
5s L = s Qp 



e «V/W 1 + j "Qp-"0n 
e qV/k B T , Sa R 1 ~ QQpQQn 

ao p 1 - a 2 0n 



1 



1 



(9) 
(10) 



The following relation connects the spin polarization 
across the depletion layer: 



OiL 



ao p (l - al n ) + Sa R (1 - a 0p a 0n ) 

1 - a 0n + 5a R ( a 0p ~ "On) 



(11) 
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Equations (JTJ, 10- will be referred to as the magnetic 
p-n junction equations. 

In the second step we wish to generalize the magnetic 
p-n junction equations to the case of several magnetic p-n 
junctions in series. Such a generalization is straightfor- 
ward in the unpolarized case, where each junction acts 
independently from the others, since the minority car- 
rier densities are fixed only by V. The inclusion of spin 
complicates the matter in the following sense. In a sin- 
gle junction 8n p , Ss p , and Ss n are the known boundary 
conditions, fully determining 5sr and ds^. Suppose we 
now connect two junctions as in the npn MBT in Fig. ^ 
Take the b-c junction to be the one in Fig. [3J Densi- 
ties n p and s p become n& e and s& e , themselves unknown, 
so that Ssr (now Ss c ) is undetermined. On the other 
hand, considering b-c to be the junction in Fig. sr 
becomes s e (s„ becomes the spin source density), and 
ni, (sl) become n& e and Sbe- These three densities are 
determined also from n p and s p , which are now denoted 
as Ufa and S{, c - This loop shows the need to obtain the 
densities inside the transistor (or a more general junc- 
tion device) self-consistently. Charge and spin are cou- 
pled both across the depletion layers [through Eqs. © 
and i)10[l] — intrajunction coupling — as well as across the 
bulk regions between two depletion layers — interjunction 
coupling. This theory is formally developed in the Ap- 
pendix. 

In the following we consider specific applications of the 
theory. Since we will deal mostly with the amplification 
mode where the excess densities in the b-c junction are 
negligible, we can get useful insights even without the 
self-consistent solutions, using only the results presented 
in this section. We refer to this as qualitative analy- 
sis. However, we support each case using a numerical 
example of a generic MBT, calculated with the full the- 
ory presented in the Appendix. The numerical model 
is a "silicon" -based MBT with the following room tem- 
perature parameters. Since the main features of the full 
theory are captured by the qualitative formulas, one can 
easily check the properties of MBTs with different param- 
eters. The parameters given below, while generic, are for 
illustration only. Unless specified otherwise, the nominal 
widths of the emitter, base, and collector are 2 /xm, 1 
/im, and 2 /im, respectively. The donor doping densities 
of the emitter and collector are Nde = 10 17 cm~ 3 and 
N dc = 
is N ab 

temperature is ni = 10 lu cm -3 . The carrier and spin re- 
laxation times are taken to be 100 ns and 10 ns, and the 
electron (hole) diffusion coefficients D n = 100 (D p = 10) 
cm 2 s _1 , all uniform throughout the sample. The carrier 
and spin diffusion lengths are L np = (Z?„ p r) 5 30 /im, 
L pn = (D pn T) - 5 « 10 fim, L sn = (D nn Ti) ' 5 « 10 fan 
and L sp « L sn . The dielectric constant is 12. We as- 
sume a spin "ohmic" contact (5s = 0) at the end of the 
collector, while at the end of the emitter an external spin 
injection gives Ss SOUICC ^ in general. 



10 cm 6 , while the acceptor density in the base 
= 10 16 cm -3 . The intrinsic carrier density at room 



IV. ELECTRICAL SPIN INJECTION 

Electrical spin injection through MBT will be studied 
in two cases: spin injection of the source spin from the 
emitter to the collector and spin injection into the col- 
lector from the equilibrium spin in the base. The spin 
injection efficiency in both cases is proportional to the 
Boltzmann factor exp(gVb e //cBT), but the physics behind 
them is rather disparate. Unless specified otherwise, we 
work in the amplification mode, where Vb e ^ ksT (for- 
ward bias) and Vb c < (reverse bias), and in the thin 
base limit, where Wb "C L n b,L s b- 



A. Source spin 

Suppose a source spin density Ss sourcc of polariza- 
tion Sa SOUICC = 8s SOUICC /N de is injected into the emit- 
ter. What is the spin response in the collector? Con- 
sider first a nonmagnetic case (ao — everywhere). 
The spin injection involves three steps, (i) The source 
spin diffusion towards b-e. At the depletion layer the 
noncquilibrium spin is Ss e — 7o l 6 e <5s source . The spin 
polarization is 5a e = 5s e /Nd e - Note that be in "fi.be 
(li,bc) means that 7* given in Eqs. ©-© are evaluated 
for the b-e (b-c) p-n junction. Since we assume that 
L se > w e , it follows that 5a e w <5a sourcc . (ii) Transfer 
of the spin into the base. From Eq. (|llfl it follows that 
a e = ab e , showing the efficiency of the spin injection by 
the majority electrons. The corresponding spin density is 
Ssbe = Sab e riQbSxp(qVbe/kBT), as follows from Eq. I|10|) 
for the forward bias case, (iii) Spin injection into the 
collector. Equation Q implies that 5s c = 7i.h c Sfce and 
so the nonequilibrium spin polarization in the collector is 
a c = "fi^bcSbe/Ndc, a result of the minority electrons spin 
pumpingisi spin in the base diffuses towards the reverse 
biased depletion layer where it is swept by the built-in 
field to the collector. Here the spin density accumulates 
as it is bottle-necked by spin diffusion and spin relax- 
ation. 

The general formula for the spin injection, combining 
the processes (i) through (iii) in a magnetic transistor 
follows from the magnetic p-n junction equations: 



5a r = Sa s 



"70,be7l,hc- 



nob eiy»e/k B T 1 _ aobaQe 

!-«0e 



(12) 



In the small injection limit nobGxp(qVb e /kBT)/Nd c is 
small (less than about 0.1); 70. be is of order one. The spin 
injection efficiency increases with increasing 71,6c. In the 
thin base limit and for a wide collector (w c 3> L sc ), for 
example, obtains ji.be ~ L sc /wb- This can be as large 
as a ten or a hundred, making a c a significant fraction 
of a e ~ a sourcc . The decrease of Wb can be achieved by 
increasing the width of the b-c depletion layer, which, in 
turn, increases with increasing |Vf, c |. 

Figure 2| illustrates the electrical spin injection of the 
source spin in our numerical model system using the full 
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FIG. 4: Calculated electron and spin density profiles (top) 
and the spin polarization (bottom) in a nonmagnetic npn 
transistor with a source spin of polarization 5a SOUICC = 0.9 in 
the emitter. The densities inside the depletion layers (shaded 
boxes) are not calculated; they are shown, with no justifi- 
cation beside guiding the eye, as straight lines connecting 
the densities at the depletion layer edges. Bias voltages are 
Vie = 0.5 volt and Vbc ~ volt. Somewhat less than 10% of 
the source spin polarization is transfered to the emitter. 



theory. The source spin of polarization <5a sourcc w 0.9 
first diffuses towards the base with a small decrease 
due to spin relaxation. The spin polarization remains 
a constant through the b-e depletion layer, resulting in 
a nonequilibrium spin density in the base. The spin po- 
larization remains steady in the base, while both n and 
s decrease to their equilibrium values in going towards 
the b-c depletion layer. Right before the depletion layer 
the spin polarization sharply drops, to get equal with 
a c w 0.075. Such sharp drops are characteristic of the 
spin pumping by the minority carriers^! 



B. Equilibrium spin 

Is there a way to accumulate spin in an MBT without 
first injecting a source spin into the structure? The an- 
swer is positive. In fact, there are two different nonequi- 
librium spin densities accumulating as a result of the car- 
rier transport through the magnetic base. The first re- 
sults from the spin extraction, acting in the emitter, the 
second from the intrinsic spin injection, effective in the 
collector. 

The extracted spin accumulates in a way similar to the 
magnetic diode4i^£ The extracted spin density is small, 
on the order of the excess minority carrier densities. It 
is necessary that the base has an equilibrium spin po- 
larization. The emitter can, but need not be, magnetic. 
Following Eq. (JTJ we get 
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FIG. 5: Spin extraction and intrinsic spin injection. Calcu- 
lated electron and spin density profiles (top) and the spin 
polarization (bottom) in an npn transistor with a magnetic 
base and nonmagnetic emitter and collector. In the emitter 
region spin density s < (extraction), here plotted as posi- 
tive (s — > — s); the spin density is not shown in the b-e de- 
pletion layer. The equilibrium spin polarization in the base is 
Qfob = 0.762, corresponding to the conduction band splitting 
of 1 ksT. The bias voltages are Vbe = 0.5 volt and Vbc = 
volt. The intrinsic spin injection, acting under the low injec- 
tion limit, results in the spin polarization in the collector of 
Sa c « 10%. 



The result is a spin extraction from the emitter, since 
the accumulated excess spin 5s e has the opposite sign 
than the equilibrium spin in the base. The extracted 
spin polarization Sa e = 8s e /Nd e is small due to the 
generally large value of N^ e . This extracted spin den- 
sity can also be treated as the spin source which prop- 
agates to the emitter region, but the contribution to 
the collector spin is negligible, being of the order of 
[n ob exp(qV be /k B T)} 2 /N de N dc . 

Intrinsic spin injection has no analog in the magnetic 
diode. The following physical processes are at work, (i) 
Minority electron injection into the base. The base has 
a spin-split conduction band, so the electrons with the 
preferred spin will move at a faster rate, resulting in a 
nonequilibrium electron minority population, but with 
the equilibrium spin polarization. The spin density is 
then out of equilibrium, (ii) The nonequilibrium spin 
density at be acts as a spin source in the b-c junction, sim- 
ilar to the spin-polarized solar cell4& (iii) This "source" 
spin is injected into the collector, where it accumulates. 

The result of the intrinsic spin injection, again in the 
limit of the thin base, is 



5a c 



ao67o,6e7i,f>c- 



n ob e 



qV b Jk B T 



N, 



(!c 



(14) 



Ss e 



-l2,be cosh(w b / L nb )s ob exp(qV be /k B T) . (13) 



If both a source , aob t ne total spin polarization in 
the collector is given by the sum of Eqs. (|12|) and 114|l . 
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FIG. 6: Charge currents in an npn MBT. The electron emit- 
ter j" and collector j™ currents are determined by the excess 
electron densities Sribe and Snbc in the base. Similarly, the 
hole emitter j% (collector j^) current is determined by 8p e 
(Sp c ). The base current jb — j e — j c is formed by the elec- 
trons recombining with holes and by the holes that both 
recombine with electrons and enter the base from the base 
electrode Shading on the arrows of the electron currents 
indicate that the current is spin-polarized. 



Remarkably, for otQ e = 0, the equilibrium spin polariza- 
tion ctQb in Eq. (|14f> plays the role of Sa e in Eq. Q14[l. 
The equilibrium spin polarization behaves, in MBTs, as 
a nonequilibrium source spin! This follows from the spin- 
selective electrical injection across the b-e depletion layer. 

Spin extraction and intrinsic electrical spin injection 
through an MBT are illustrated in Fig. [5] using the full 
theory. The equilibrium polarization spin in the base is 
Q'ofa = 0.762. The electrical transport through the base 
leads to a spin extraction from the emitter, with the ex- 
tracted spin polarization a e « —0.001, small due to the 
large value of Nd e . The spin polarization jumps to its 
equilibrium value in the base, increasing sharply (see the 
discussion to Fig.QJ to a^ b + 8a c right before reaching the 
second depletion layer. The injected spin polarization is 
5a c pa 10% , relatively large due to the small value of Nd c 
and the large ratio L sc /w b [see Eq. 

We expect that both the direct injection of the source 
spin as well as the spin extraction and the intrinsic spin 
injection become more efficient in the limits of large car- 
rier injection (large biases), where our theory does not 
apply. This expectation is based on the results of numer- 
ical calculations^ 7 of spin injection in magnetic diodes. 



V. ELECTRICAL CHARACTERISTICS 

The electrical properties of MBTs are determined by 
both the charge and the spin of the current carriers. 
There are two ways spin affects the electrical currents: 
through the dependence of the equilibrium minority elec- 
tron and hole densities on the equilibrium spin polariza- 
tion, and through the spin-charge coupling resulting from 
the presence of a nonequilibrium spin. We first introduce 
the formalism for calculating electrical currents in bipolar 
transistors and then analyze the two ways in detail. We 
conclude with a discussion of the spin current through 
MBTs. 



The scheme and the sign convention for the currents 
is shown in Fig. [5] (see also Fig. ^ for the description 
of symbols labeling the regions). Below we summarize 
the expressions for the currents from the Shockley the- 
ory of bipolar transistors We write the expressions 
in a rather general form which turns out to be applicable 
also to MBTs (this follows from the generalized Shockley 
theory of magnetic diodes^i). Let us define the genera- 
tion current of a carrier c (electron or hole) in region r 
(emitter, base, or collector) as 



igr 



qD cr ( W r 

— — c 0r coth I — 

±Jr r V J-Jri 



(15) 



The electron charge current density in the emitter is 



Je 



Jgb 



Sn be 1 Sn bc 

n ob cosh(w b /L nb ) n ob 



(16) 



The first term in Eq. 1)16(1 represents the diffusion of elec- 
trons in the base at the b-e junction. The second term 
describes a competing diffusion from the excess minority 
electrons at the b-e junction. Through the current conti- 
nuity, this base diffusion current continues to the emitter 
to become J™. Similar expression holds for the electron 
current in the collector: 
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be 
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n ob cosh(w b /L nb ) n ob 



(17) 



All the carrier densities appearing in the expressions for 
the currents can be calculated from the theory in Sec. lIIII 
Holes contribute to the currents through the diffusion 
of their excess minority populations Sp e and Sp c . The 
expression for the densities is similar to Eq. Q with no 
spin. For example, 



Sp e = POe (t 

The two hole currents arc 
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(18) 
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(19) 
(20) 



POe 
Spcb 
POc ' 

The total emitter current is j e = and similarly the 

total collector current is j c = j™ + j%. Finally, the base 
current is j b — j e ~jc- The task of computing the currents 
through an MBT is reduced to the computation of the 
excess electron and hole densities at the two depletion 
layers. 

The current amplification coefficient (gain) (3 is defined 

as 







J_c_ 

Jb ' 



(21) 



If (3 is large, small changes in j b lead to large variations 
in j c , allowing signal amplification. The gain is often 
written as^i 



/3=l/K + 7 '). 



(22) 



Here a' T measures the (in)emciency of the electron- 
hole recombination in the base, and 7' describes the 
(in)efhciency of the emitter electron injection into the 
base. The usual base transport factor— is defined as 
dT = 1/(1 + ot' T ) and the emitter efficiency factor as 

7 = 1/(1+7')- \ 

After substituting for the currents, the emitter 

(in) efficiency a' T is calculated to be 



500- 



cosh(- — ) 

L n b 



1. 



(23) 



This is the value obtained for conventional transistors 
and is also valid for MBTs. Coefficient a' T does not de- 
pend on spin since it reflects only the electron-hole re- 
combination in the base and in our model the recombi- 
nation is spin independent. However, there may be cases 
where L nb depends significantly on ao b , in which case the 
gain could be controlled by spin even through a' T . Note 
that Wb depends on the equilibrium spin through the spin 
dependence of the built-in fields. 41 Eq. i|23|) holds even 
for such cases. On the other hand, 7' depends explicitly 
on both the equilibrium and nonequilibrium spin. We 
describe this dependence by defining a new parameter^: 



i = 1'0/v, 



(24) 



where 7q is the emitter efficiency of a conventional npn 
transistor— 



7o 



N ab D pe n 2 ie L nb smh(w b /L nb ) 
N de D nb nf b L pe tanh(w e /L pe ) ' 



(25) 



where we allow for a generally different intrinsic carrier 
concentrations n ie and n ib in the emitter and base, re- 
spectively. In the next two sections we discuss the physics 
behind ?y, which we call the magnetoamplification coeffi- 
cient. We will in particular consider the thin base limit, 



where 



(w b /L nb ) 2 , 7q ~ (wb/Lnb), and j' domi- 



nates the current amplification (for example in Si). In 
such cases 

P^vh'o- (26) 

If the base transport factor is not negligible, the spin 
control efficiency diminishes. 



A. Magnetoamplification effect: influence of the 
equilibrium spin 



Consider a magnetic base. The Boltzmann statistics 

41 47 

give s 1 ' 



n ob 
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N al 



tanh(q( b /k B T) 
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Nab y/T 



(27) 
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Since j™ 6 ~ n ob [see Eq. <jl5j) ]. it follows that the base gen- 
eration current increases as ao b (that is, its magnitude) 
increases. In turn, j™ , ~ j™ b , so that both the emitter 
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FIG. 7: Calculated gain /3 of an MBT with a magnetic base 
(solid) and emitter (dashed) for our numerical model. No 
source spin is present. The equilibrium base (emitter) spin po- 
larization is aot (aoe). The calculation is done on a structure 
with the nominal base width of 1.5 /xm to keep the effective 
width wt positive for the considered range of polarizations. 



and the collector currents increase with increasing a^ b . 
The equilibrium spin controls the charge currents flow- 
ing in MBT, leading to a magnetoresistance effect. Spin 
unpolarized holes too contribute to the spin control of 
the currents, as shown by Lebedeva and Kuivalainen^ 
for a pnp MBT. If the emitter is magnetic, the minority 
hole density is 



1 



PCX: 



N dc y/T 



(28) 



analogously for the collector (in the amplification mode 
the hole density in the collector is negligible and does not 
affect the current properties). The hole emitter current 
increases with increasing ctQ e . 

The magnetoamplification coefficient amounts to the 
simple expression 



V = 



Hie 



(29) 



and in the thin base limit, Eq. I|26|l is thus described by 
the gain 




(30) 



'Ob 



The gain of MBTs can thus be controlled by controlling 
the equilibrium magnetization (for example, by changing 
the external magnetic field) of the emitter or the base. 
The collector magnetization plays no role. On the other 
hand, the two equilibrium polarizations ao e and ao b act 
against each other: The gain increases (decreases) with 
a greater spin splitting in the base (emitter). This is 
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because the emitter (in)efficiency 7' increases (decreases) 
if there are relatively more holes (electrons) present in j e . 
If the spin polarization is uniform across the b-e junction, 
the gain is spin independent. 

The opposite role of the equilibrium magnetizations in 
the base and in the emitter is shown in Fig. which il- 
lustrates the behavior of (3 with respect to the changes of 
ctQb and aoe separately. The calculation is done using the 
full theory, not the approximate formulas above. How- 
ever, the approximation describes the calculation very 
well, showing that spin dependent effects for example on 
the effective widths w, which are accounted for in the full 
theory, play minor role in our example. 



B. Giant magnetoampliflcation effect: spin-charge 
coupling in MBT 

A nontrivial realization of the Silsbee- Johnson spin- 
charge coupling representing the physics of the prox- 
imity of an equilibrium and nonequilibrium spin in MBTs 
is what we call here the giant magnetoampliflcation ef- 
fect (GMA), in analogy with giant magnetoresistance 
(GMR) effect in metallic multilayersi2 For GMA it is 
necessary that there be a nonequilibrium spin polariza- 
tion in the emitter (arising from a source spin) and an 
equilibrium spin either in the base or in the emitter (or 
both) . The physics is illustrated in Fig. [21 The charge 
current through MBTs depends on the relative orienta- 
tion of the source and the equilibrium spins, because of 
the spin-dependent barrier in the b-e junction. 

In the presence of a nonequilibrium spin density 5a e , 
the magnetoampliflcation coefficient becomes 



V = 



1 -<*0e 



[l + 5a e (a Qb - aoe)/(l - a 0e ) 2 ] ■ (31) 



If only the base is magnetic, the gain in the thin base 
limit is 



= 



11 + Sa e a b 



at, 



(32) 



01, 



The spin-charge coupling is described by the product 
Sa e aob, similar to implications of the spin-voltaic effect 
in magnetic p-n junctions4i*2LSi Let /3 max and /3 m j n are 
the gains for the configuration of the source and equilib- 
rium spins (parallel or antiparallel) that yield the max- 
imum and minimum gain, respectively. For a magnetic 
base (emitter) the maximum is achieved at parallel (an- 
tiparallel) orientation and the minimum at antiparallel 
(parallel) orientation of the source and equilibrium spins, 
respectively. We define the GMA coefficient as 



GMA 



max Mmm 



0B 



(33) 



in analogy with a similar expression (involving resistivi- 
ties) defining the GMR coefficient. 



For the magnetic base 



GMA 



2\Sa e aob\ 
1 - |<5a e aob| 



(34) 



If, for example, Sa e = a Qb = 0.5, GMA=67%. The 
analogy with GMR is clear: there is a large magnetore- 
sistance effect (greater than 10%), which is most pro- 
nounced when the relative orientation of two spin po- 
larizations changes from parallel to antiparallel. If, on 
the other hand, the emitter is magnetic, the effect is op- 
posite: the parallel spin orientation decreases the gain, 
due to the decrease in the the emitter injection efficiency. 
The GMA coefficient is 



GMA 



2|(5a e aoe 



1 - ai 



\Sa e a 0e 



(35) 



The GMA coefficient vanishes if a^ e = a^b- To decide 
on whether to use a magnetic base or a magnetic emit- 
ter one needs to take into account that a magnetic base 
will have a smaller L s b- If L s b < Wb, a magnetic emitter 
would be instead preferable. Note that the GMA coef- 
ficient is directly proportional to the magnitude of Sa e , 
and so it can be used to measure the nonequilibrium spin 
polarization, as in the case of magnetic diodes^i 

Figure il illustrates GMA on our MBT example with a 
magnetic base and source spin polarization a sourcc = 0.9 
(as in Fig. EJ) . The solid line represents the calculation 
of f3(ctob) using the full theory, while the dashed line is 
the approximation Eq. (|32() valid for the thin base tran- 
sistors. The approximative formula works very well. The 
asymmetric curve demonstrates the GMA effect. When 
the equilibrium and the source spin are antiparallel, [3 is 
small; when they are parallel, j3 is large. The magnetic 
control of the charge current and its amplification is thus 
predicted to be rather effective. 



C. Spin-polarized current 

Thus far we have studied the spin polarization of the 
electron density as a measure of the spin injection effi- 
ciency. This spin polarization is typically detected by 
optical experiments. Another spin polarization, that of 
the charge current, is more convenient for theory and is 
invariably used in spin injection models to assess the spin 
injection efficiency (see the 7-technique of Rashba 52-53 ) 
and to establish the boundary conditions for spin at the 
interfaces^ The current spin polarization ay is the ratio 
of the spin current j s , which is the difference between 
the charge currents formed by the spin up and spin down 



electrons, and the total charge current j: 



■ js/j- The 



current spin polarization is much less intuitive than the 
density spin polarization. There are several reasons for 
that. First, etj involves not only the carrier and spin den- 
sities, but also the drift or diffusion velocities which can 
be spin dependent. Next, unlike a, the magnitude of ctj 
is not restricted to the interval (0, 1). The magnitude of 




FIG. 8: Giant magnetoamplification effect. Calculated gain /3 
of our MBT example with a magnetic base and source emitter 
spin, as a function of aob for a fixed Q sourco — 0.9. The 
solid curve is the calculation using the full theory, the dashed 
curve is the contribution of the emitter efficiency only, /3 ps 
77/7', where the magnetoamplification coefficient 77 is given by 
Eq. 13111 . The calculations are done on a structure with the 
base long 1.5 /im, as in Fig. Q 



dj is not even bound. Indeed, the spin current can be 
finite even if there is no charge current, making aj infi- 
nite. Also the signs of aj and a can be different. Finally, 
unlike the charge current, the spin current need not be 
uniform. Because of spin relaxation, j s (and also ay) is 
not conserved. For the above reasons, unless the relation 
between a and aj is obvious or is explicitly derived, aj 
is not indicative of the spin injection efficiency. In par- 
ticular in inhomogeneous (or hybrid) semiconductors at 
degenerate doping densities or at large biases, both spin 
diffusion and spin drift are relevant and one needs to em- 
ploy the Poisson equation to solve the transport problem 
self-consistently to obtain a dependence between a and 
aj. Many of the experimental spin injection results are 
likely to fall in this category, making realistic theoretical 
modeling difficult. 

Fortunately, in the low injection limit there is a sim- 
ple relation between the spin current and the spin, so 
the knowledge of aj together with the knowledge of the 
charge current j suffice to obtain a. The spin current 
density at point c (see Fig. ^| is readily obtained from 
Tab. II in Ref. El 
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(36) 



Thus, Sa c = Ssc/Ndc is directly proportional to j sc = 
jaj. As is shown below, aj is usually comparable to 
aQ or 8a, largely independent on the biases. The spin 
injection efficiency is then determined by j c , which, in 
turn, depends exponentially on V& e . 

We adopt the same sign convention for the spin cur- 
rents as for the charge currents, see Fig. EI A straightfor- 
ward application of the magnetic p-n junction equations 



FIG. 9: Calculated spin current polarization profile in the 
magnetic npn transistor, calculating by the full theory. The 
source spin at the emitter is fixed at <5a sourcc = 0.9, while the 
equilibrium base spin polarization changes sign at the magni- 
tude of Qob = 0.762, corresponding to the base spin splitting 
of the order of the thermal energy. The biases are V& e = 0.5 
volt and Vbe = volt. The spin current remains uniform 
across the depletion layer, based on the assumption of neg- 
ligible spin relaxation there. The spin current is normalized 
by the emitter and the collector currents in the emitter and 
the collector regions, and by the base electron current in the 
base. The horizontal line is for no source spin. The spin cur- 
rent stays uniform throughout the sample at the value given 
by aob, implying no significant spin relaxation, but not nec- 
essarily spin injection. 



fSec. IHIjl gives 
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What is striking (although not so surprising) is that aj c 
in the cases of practical biases |Vj, e | ksT is bias inde- 
pendent. Indeed, for a forward bias Vb e , that is, in the 
amplification mode, 



a b(l - cvqJ + 5a e (l - a bOto e ) 



1 



+ 6a e (a b - «oe) 



(38) 



This is in sharp contrast, for example, to Eq. (|12|) which 
displays the exponential increase of the spin injection ef- 
ficiency with Vb e ■ In the limit of no source spin, Eq. I|38|l 
reduces to aj C — aob, the main result of Flatte et ali2& 
Similarly, if the transistor has no equilibrium spin, the 
spin current polarization is a JC = 8a e . If the emit- 
ter/base bias is reverse, Vbe < 0, spin injection is prac- 
tically noncxisting. Yet, there is a large current spin 
polarization, <Xj C = arjfc, independent of the source spin, 
confirming our claim that large aj alone does not imply 
efficient spin injection. 

Our analysis is illustrated in Fig. |5J We consider the 
amplification mode, but the results vary little with the 
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applied bias, as explained above. Figure|5]describes three 
cases. One with the source spin and the equilibrium spin 
pointing in the same direction, one where the two spins 
are antiparallel, and one where only the equilibrium spin 
is present. In all the cases otj c (which is the value at 
x ~ 3 /im) is smaller than 1, but rather considerable, 
consistent with Eq. (|38fl . Note that ay in the emitter 
is much greater than 1 for the case of the source spin, 
due to the spin diffusion current being greater than the 
majority electron drift current there. The value of ay 
decreases at the emitter/base junction, where it becomes 
a constant, signifying small spin relaxation in both the 
base and in the emitter. Finally, the case of no source 
spin polarization {8a e <C aof>) shows a flat ay. This case, 
viewed from the density spin polarization perspective, is 
shown in Fig. [SJ The nonequilibrium spin which accumu- 
lates in the base is injected to the emitter, and extracted 
from the collector. The injection/extraction efficiency in- 
creases exponentially with Vf, e . For small biases, even if 
ay is equal to aob ~ 0.762, the spin injection is negligi- 
ble. Interestingly, even if the spin current polarization is 
positive in the emitter, the spin density polarization is 
negative (spin extraction, cf. Fig.[SJ). 



VI. CONCLUSIONS 

We have developed an analytic theory for the spin- 
polarized transport in magnetic bipolar transistors in the 
low injection regime. We have shown that the transistor 
displays a number of novel phenomena, not observed ei- 
ther in the conventional spin unpolarized bipolar transis- 
tor or in the magnetic diode. One such effect is the intrin- 
sic electrical spin injection, which is a spin injection from 
a magnetic base into a nonmagnetic collector, without 
the presence of a source spin in the emitter. The spin 
injection efficiency depends exponentially on the emit- 
ter/base bias. Other effects are related to the gain which 
can be influenced in two ways. First, by the equilibrium 
spin, either in the emitter or in the base, modifying the 
generation current. Second, by the spin-charge coupling, 
modifying directly the electron injection from the emit- 
ter. 

Bipolar junction transistors are used for ultra high 
speed logic applications thanks to the fast carrier trans- 
port. Spin can bring additional functionalities. The mag- 
netoamplification effects can be used to study spin sig- 
nals (time varying spin polarizations) by detecting charge 
currents. Indeed, all the currents in an MBT in the am- 
plification mode depend primarily on Sribe, which is con- 
trolled by both the equilibrium and nonequilibrium spin. 
The spin-to-current conversion can thus be observed by 
measuring the collector current, or directly the GMA co- 
efficient Eq. (|33[) . On the other hand, the amount of 
the injected spin in the collector depends exponentially 
on Vbe in the same manner as the charge current de- 
pends on it. Changes in, for example, the base current, 
can thus lead to changes in the spin polarization 5a c . 



As a result, current signals can be detected by observ- 
ing the accumulated spin polarization. Perhaps the most 
attractive feature of the magnetoamplification effects is 
that the spin splitting or the source spin polarization are 
not built-in device properties but can change on demand, 
during the transistor operation, by magnetic field. This 
is why an MBT is an example of a variable heterostruc- 
ture transistor. 

Another use of MBTs may be in the electrically in- 
duced magnetization switching, similarly to what has 
been observed as light-induced ferromagnetismi2ii4 or 
ferromagnetism induced by the gate voltage of field-effect 
transistors If the base is a ferromagnetic semicon- 
ductor, the equilibrium magnetization depends on the 
density of free carriers. This density can be, in turn, 
controlled by \%e- While the scenario of the electri- 
cally induced ferromagnetism in a nondegenerate MBT 
is probably not realistic, at higher doping and current 
injection levels (where our theory no longer applies) this 
effect could be observable. 

We believe that the phenomena we propose to study 
are robust and should be observed. For the spin-source 
spin injection one does not need a magnetic semiconduc- 
tor in the structure. The source spin can be generated in 
the emitter either optically or electrically, and similarly 
the spin injected into the collector can be observed by 
detecting the polarization of electroluminescence^^ For 
the phenomena related to the equilibrium magnetism, the 
rule of thumb is that the spin splitting should be compa- 
rable to the thermal energy in order to create noticeable 
spin polarizations. While to see GMA at room tempera- 
ture may be difficult at present, at smaller temperatures 
(starting from 50 K where the shallow donors and accep- 
tors start to ionize) the effects could be detected at the 
spin splitting levels of ~ 1 meV. 
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APPENDIX A: AN ARRAY OF MAGNETIC P-N 
JUNCTIONS 

We introduce a formalism for evaluating the carrier 
and spin densities in an array of magnetic p-n junctions. 
An array of two junctions forms an MBT, while three 
junctions would form a magnetic thyristor. The bound- 
ary conditions are applied for the densities; it is straight- 
forward to adapt the method to have boundary condi- 
tions specified by the spin currents. All the junctions 
need to be considered simultaneously since there is both 
intra- and inter- junction charge and spin coupling. The 
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gion as a result of the carrier injection: 
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FIG. 10: Notation scheme for an array of p-n junctions. The 
figure illustrates an npnp (thyristor) structure. The depletion 
layers are shaded. The known quantities are shown above the 
junctions, the unknown below. Each junction i is character- 
ized by the scalar and by the triplet of vectors v°, Ci, 
and Di, which are determined by the doping densities and 
the equilibrium carrier and spin densities of the regions adja- 
cent to the junction, and by the applied voltage Vi across the 
junction. The nonequilibrium spin density uo (note that the 
symbol here denotes the region, not an equilibrium) and the 
charge/spin density vector V4 would be the boundary condi- 
tions here. The densities u% and v< at the depletion regions 
are to be obtained self-consistently. These densities are all 
what is needed to calculate the charge and spin currents. 



intrajunction coupling arises from the generalized Shock- 
ley conditions— of the uniformity of the spin-resolved 
chemical potentials and of the continuity of the charge 
and spin currents, across a junction's depletion layer. 
The interjunction coupling arises from the carrier and 
spin diffusion in the bulk regions between the depletion 
layers: the carrier and spin density at one end influences 
the current at the other end, and vice versa. The two 
couplings lead to a set of linear algebraic equations for 
the densities. 

The array we consider is shown in Fig. 1101 Each junc- 
tion is given a number i starting from one. The junctions 
are either of the p-n (p left and n right), or of the n-p (n 
left and p right) type. In the following the indexes n and 
p relate to the n and p regions adjacent to the junction in 
question. Let us introduce the notation using a generic 
junction, as in Sec. IIIII The bias V across the junction is 
positive for the forward and negative for the reverse di- 
rection of the charge current. The electron density is n, 
with index zero (no) reserved for the equilibrium value. 
The nonequilibrium (excess) part of the density is de- 
noted as 5n = n — rig. Similarly for the electron spin 
density s and for the spin polarization a = s/n. The 
doping densities are Nd for the donors in the n region 
and N a for the acceptors in the p region. 

We make the complex notation more compact by in- 
troducing some unifying symbols. We first denote by u 
the scalar characterizing the nonequilibrium spin density 
due to the carrier extraction: 

u° = - 72 cosh(w p / L np )s 0p (f v / kBT - l) . (Al) 

Here fcgT is the thermal energy, with T denoting tem- 
perature; q is the proton charge. Parameters 7 are intro- 
duced in Eqs. To properly account for the spin 
current in the n region, one needs to apply the rescaling 
in Eq. (jSJ). We also introduce vector v° which specifies 



(A2) 

The dimensionless vector C characterizes intrajunction 
coupling: 

c = K)p(72 ~ 7i), 7i] > (A3) 
while another dimensionless vector D, given by 

nop • 



D 



e qV/k B T 



N d 1 - a% n 



[a Qp ~ a 0n , 1 - a 0p a 0n } 



(A4) 



characterizes interjunction coupling. The quantities u , 
v°, C, and D are presumed to be known. They are the 
input parameters. 

The unknown quantities are the nonequilibrium carrier 
and spin densities at the depletion layers. In the low- 
bias regime considered here, the electron density in the 
n regions is fixed: n n « Nd, and only Ss n needs to be 
calculated. In the p regions the excess electron (spin) 
density Sn p (Ss p ) is unknown. We denote by scalar u the 
nonequilibrium spin density 8s n in the adjacent n region: 



5s n . 



(A5) 



Vector v will describe both the nonequilibrium electron 
and spin density in the p region: 



v = [Sn p , 8s p ] 



(A6) 



The boundary conditions are given by the corresponding 
u or v (depending on whether the n or the p region is 
the contact region) at the beginning and the end of the 
array. 

Using the above notation, the magnetic p-n junction 
equations Eqs. (JIJ, (J3J, and ifTUI) are greatly simplified 
and can be adapted to solve the array problem. Indeed, 
for junction i the p-n junction equations can be written 
as 

Ui = m° + 7o,»u J ±i + C, • v iT i, (A7) 
v, = Vl °+D lUl , (A8) 

where the upper (lower) sign is for the p-n (n-p) directed 
junction. For the npn MBT described in the main text, 
the equations take the form 



wi = u 1 +j . 1 u + C 1 -v 2 , 
vi = v° + Dim, 



for junction 1, and 

U 2 = ^270,2^3 + C 2 ■ Vi, 



v 2 



D 2 it 2 , 



(A9) 
(A10) 

(All) 
(A12) 



for junction 2. Analogous equations can be written for 
holes. The solution to Eqs. (|X9|) - (|A12|) is 



"2 



(C 2 ■ Di) (70,11*0 



C 2 -v° 



-70.2M3, 
(A13) 
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where we have neglected the terms of order 
[n 0p exp(V)/N d ] 2 , small in the low bias regime. The 
formulas for u\, vi, and V2 can be obtained directly by 
substituting Eq. (|A13j) back to Eqs. I|A9|I through l|A12|l . 

Equation l|A13() describes spin injection through an 
MBT, since U2 is the nonequilibrium spin in the col- 
lector at the depletion layer with the base. The first 
term on the right-hand side (RHS) of Eq. (|A13|I repre- 
sents the transfer of source spin uq from the emitter to 
the collector. Indeed, for a nonmagnetic transistor (the 
equilibrium spin polarizations are zero) the transferred 
source spin is 113 — 7o,i7i,2«Op exp(Vi/A:BT)uo- Here 
70 describes the transfer of the source spin through the 
emitter — a majority carrier spin injection. Once the spin 
is in the base, it becomes the spin of the minority carriers 
[hence the minority density factor nob exp(Vi/fcsT)], dif- 
fusing towards the base/collector depletion layer 2. The 
built-in electric field in this layer sweeps the spin into 
the collector, where it becomes the spin of the major- 
ity carriers again, by the process of the minority-carrier 
spin pumping4i^& Can the injected spin polarization in 
the collector be greater than the source spin polariza- 
tion? The answer is negative in the low-injection regime. 
It would be tempting to let the spin diffusion length in 
the collector to increase to large values to get a greater 
pumped spin. But that would increase the importance of 
the electric field in the n-regions and the theory (which is 
based on the charge and spin diffusion and not on the spin 



drift) would cease to be valid. However, the spin density 
in the collector can be greater than that in the base, 
demonstrating that the spatial decay of the nonequilib- 
rium spin (spin accumulation) is not, in general, mono- 
tonically decreasing, in line of what was demonstrated 
nonmagnetic p-n junctions 4 *?- as well as inhomogeneous 
unipolar semiconductors^* The same term also describes 
the transfer of the nonequilibrium spin 11® accumulated as 
a result of the electron injection into the magnetic base. 
The second term of Eq. I|A13|) , which is independent of 
the source spin, results from the (intrinsic) spin pump- 
ing by the minority channel of the nonequilibrium spin 
generated in the base by the forward current through 
junction 1. This term vanishes if the base is nonmag- 
netic (aoh = 0). The third term of Eq. i|A13(l repre- 
sents the spin extraction due to the magnetic base. This 
term is controlled by V2, the bias at the base/collector 
junction. For the reverse bias, used in the amplification 
mode, this term can be neglected. Finally, the last term 
of Eq. I|A13I) , describes the diffusion of the source spin in 
the collector. 

The knowledge of the carrier and spin densities at the 
depletion layers allows us to calculate the charge and spin 
currents in the systems, as well as the density spatial 
profiles throughout the bulk regions, using the formulas 
for magnetic p-n junctions in Tab. II of Ref. This is 
done for our numerical MBT model presented in Sees. II VI 
and E| 
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